o 

(N 



Quantum memory for light via stimulated off-resonant Raman process: 
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We consider a quantum memory scheme based on the conversion of a signal pulse into a long-lived 
spin coherence via stimulated off-resonant Raman process. For a storing medium consisting of alkali 
atoms, we have calculated the Autler-Townes resonance structure created by a strong control field. 
By taking into account the upper hyperfine states of the D\ optical transition, we show important 
deviations from the predictions of the usual three-level A-scheme approximation and we demonstrate 
an enhancement of the process for particular detunings of the control. We estimate the memory 
efficiency one can obtain using this configuration. 

PACS numbers: 42.50.Gy, 42.50.Ct, 32.80.Qk, 03.67.-a 
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I. INTRODUCTION 

Long-lived and highly efficient quantum memories for 
fight are a crucial tool for quantum information process- 
ing, quantum computing, long-distance secure commu- 
nication and scalable networks [I] [5] . Various physical 
systems are intensively investigated and many significant 
advances have been recently reported [3j i4j . In particu- 
lar, large ensembles of identical atoms, as gases at room 
temperature [5'~'5^ or ultracold samples [TUin^ , have been 
successfully used for demonstrating the storage and re- 
trieval of quantum states in different regimes. 

In the theoretical investigations of such quantum mem- 
ories, the description of the light-atoms interface is usu- 
ally based on a three-level A configuration, with two 
ground states and one excited state [I3til5j . However, the 
hyperfine interaction in the upper states of alkali atoms 
is not strong enough for the system to be perfectly de- 
scribed by this approximation. The multilevel structure 
of alkali atoms has thus to be taken into account for a cor- 
rect description of the quantum memory effect [161 I17j . 
In this paper, we show for the Di optical transition how 
an additional level significantly modifies the stimulated 
Raman process and we finally evaluate the efficiency of 
pulse storage and retrieval in this configuration. 

The paper is organized as follows. Section |Tl] first in- 
troduces the energy diagram and the assumptions on 
which our work is based. We then present the theoreti- 
cal model and provide a general expression of the atomic 
sample susceptibility. In Section |III[ we study the sus- 
ceptibility behavior for different frequency detunings of 
the control and probe modes and we explain how the 
multilevel structure significantly modifies the quantum 
light transport driven by coherent scattering mechanism 
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in the conditions of electromagnetically-induced trans- 
parency or stimulated Raman process. In section |IVJ we 
present calculations for the transport of a signal pulse in 
an optically dense sample. The signal pulse is tuned in 
the spectral region where the Autler-Townes resonance 
would be maximal in its amplitude and the stimulated 
Raman scattering would be most effective. The achiev- 
able efficiency of a quantum memory protocol performed 
in this configuration is finally discussed, in both cases of 
forward and backward retrieval. 



II. AUTLER-TOWNES EFFECT IN THE 
Di-LINE OF ALKALI ATOM 

In this section, we present the energy configuration 
under study and we detail the theoretical model used to 
fully describe this configuration. General expression of 
the susceptibility of the medium is finally given. 



A. Basic assumptions 

In this paper we consider the case of alkali atoms and 
focus in particular on the Di optical transition. Figure [T] 
gives the energy diagram and figure [2] shows the excita- 
tion geometry. The atoms populate the upper hyperfine 
sublevel of their ground state with maximal spin projec- 
tion {F^ — I + 1/2, M = F_|_} (where / is atomic nuclear 
spin) and we denote this state as \m). 

The control field has a right-handed circular polar- 
ization (cr+). In this configuration, there is no inter- 
action of the control field with the populated sublevel 
and only the presence of the probe mode in the or- 
thogonal left-handed circular polarization (cr_ ) opens the 
excitation channel. The probe and control modes co- 
herently couple the populated state |m) with the Zee- 
man sublevel {F^, M = F+ ~ 2} in the ground state. 
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FIG. 2: (Color online) Schematic diagram showing the exci- 
tation geometry for the transitions specified in figure [T] 

than either the natural decay rate 7, the Rabi frequency 
of the control mode or the detuning A. We will show 
the limit of this widely-used approximation via the full 
calculation of the AT structure for the Di-line of ^'^'^Cs 
atom, which has the largest hyperfine splitting among 
alkah atoms flHF = 1168 MHz = 2567. 



B. Calculation approach 



FIG. 1; (Color online) Energy levels and excitation channels 
considered for the Di line of alkali atom. The atoms populate 
the upper hyperfine sublevel of their ground state. The sys- 
tem is dressed by a strong control field with a right-handed 
circular polarization ((t+) and a frequency detuning A. The 
system is probed by a weak mode with left-handed circular po- 
larization {(J-) and frequency detuning A. The three Autler- 
Townes resonances (AT-triplet) are shown by the red dashed 
lines. 



which we denote as \m'), and with two Zeeman sublevels 
in the excited state {F!_ = I - 1/2, M' = F+ - l} and 
{F^ =7+1/2, A/' = F+ — 1}, which we denote respec- 
tively \n) and \n'). 

Such a configuration leads to the Autler-Townes (AT) 
effect [IHl [H] ■ As shown in figure fl] the energy levels of 
the excited atomic states \n) and \n) are modified by the 
presence of the control field. If the control field is tuned 
precisely in resonance with the non-disturbed transitions, 
either \n) or |n'), the interaction with the control mode 
splits the originally degenerate system of the atomic and 
field oscillators into two quasi-energy sublevels. How- 
ever, the interaction of the control with the other upper 
hyperfine atomic state is not negligible and the presence 
of the control mode affects on the locations of both the 
hyperfine resonances. For an arbitrary detuning of the 
control, we refer to the joint atom-field quasi-energies as 
the AT-triplet (instead of the usual AT doublet) reso- 
nance structure. The location of these three resonances 
for a particular detuning of the control mode is shown by 
the dashed red lines in figure [T] 

Let us note finally that the A-schcmc approximation is 
the asymptotic limit of our approach: in this limit, the 
hyperfine interaction is considered to be much stronger 



We now turn to our theoretical model. The AT reso- 
nance structure can be observed by passing through the 
sample a probe pulse that will excite the atoms from the 
populated state \m). The propagation of the probe mode 
through the medium, in the geometry shown in figure [21 
is described by the standard macroscopic Maxwell equa- 
tion, where the dielectric susceptibility of the medium, 
being linear in response to the probe mode, accumulates 
all the orders of nonlinearities for the control mode. Due 
to the linearity for the probe field, we can consider an 
arbitrary set of probe field frequencies forming the signal 
pulse and we follow the temporal and spatial dynamics of 
this pulse in a finite sample. The dynamics of the signal 
pulse and control mode come from the characteristics of 
the stimulated Raman process, as shown in [111 [20] • The 
absorption and dispersion parts of the AT-resonances in 
the sample susceptibility result from the combined action 
of this process on the signal pulse transport. 

The positive frequency component £^^' (r, t) of the sig- 
nal pulse can be described by a wavepacket propagating 
along the direction z, with a carrier frequency Co and 
wavenumber k = Co /c, such as: 



£'^+\r,t)^e{z\t)e 
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(2.1) 



For the sake of simplicity the transverse profile of the 
wavepacket is not taken into account, and we thus ig- 
nore any diffraction and transverse inhomogeneity of the 
atomic sample during the propagation process. However 
let us point out that the analysis could be generalized for 
an inhomogeneous system as well. 

We then make use of the Fourier expansion of the slow- 
varying amplitude e(z; t) of the probe field: 
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(2.2) 



The Fourier component e(z; il) obeys the following 
Maxwell equation: 
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where x{z;^) is the Fourier component of the sample sus- 
ceptibility. The spectral dependence of the sample sus- 
ceptibility is shifted by the carrier frequency such that: 



X{z:t~t')^e^^^'-''h{z:t-t') 



(2.4) 



where xi^'i^ ~ t') is the kernel of the susceptibility op- 
erator standardly defined in time representation and de- 
termining response of the polarization wave at time t on 
the driving signal wave (2.1) considered at retarded time 
t' < t. 

The susceptibility of the medium is then given by the 
following expansion: 



xiz;n) 



E E ^ (de):,™ (de) 



ni^7i,7i' 712— n.n' 



(2^ 



-"-o(2)/o(p) 



xGl;-)(pi,p, + ;ifc;a(cS + 0) + |-) (2.5) 

We have used the following notations: the transition ma- 
trix elements of the dipole operator d are projected onto 
the polarization vector of the probe mode e, ri,o(z) is the 
density distribution of atoms, and /o(p) is their momen- 
tum distribution, normalized such as 
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■/o(p) = 1. 



(2.6) 



The spectral behavior of the susceptibility (2.5) is de- 



termined by the spectral properties of the contributing 
Green's functions Gnin2 ■, describing the dynamics of the 
upper atomic states dressed by the control mode. The 
derivation details and precise definitions of the Green's 
function, matrix elements etc. are given in Appendix [K\ 
We point out here that these functions for all possible 
, 7i' and T7-2 — n,n' built the 2x2 matrix so that 
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the expansion (2.51 consists of four terms. Two diago- 



nal terms asymptotically (for highly resolved hyperfine 
structure) reproduce the A-scheme results for either up- 
per state \n) or |n'). The two others terms, which are 
interference terms, significantly modify the predictions 
of the A-scheme approximation, as we will show in the 
following. 



III. SPECTRAL BEHAVIOR OF THE SAMPLE 
SUSCEPTIBILITY 

In the case of monochromatic excitation, it is conve- 
nient to analyze the Fourier component of the dielectric 



susceptibility as a function of either its carrier frequency 
u) or its detuning A = uj — ujq from the atomic resonance 
frequency ujq. In this section we present our calculation 
of the sample susceptibility x = x(^) = x'{^) +*x"(^)i 
which were obtained for an homogeneous medium con- 
sisting of ^•^^Cs atoms. We assume a typical configura- 
tion of magneto optical trap (MOT) such that the effects 
of atomic motion on the time scale of a few microseconds 
can be disregarded. 



In the following, the susceptibility given by Eq.(2.5) 
will be scaled in units of no A , where uq is the density of 
atoms and A = A/27r . In typical MOT conditions, the 
atomic gas is dilute and tiqA <^ 1 but the sample can 
be optically thick if hqX L ^ 1, where L is the sample 
length. It is very important to underline that the tran- 
sition matrix elements Vnm' and Vn'm', contributing to 
the Green's functions given by Eqs. (A4), (A5|, cannot 



be considered as independent parameters. They are ex- 
pressed by one reduced matrix element multiplied by two 
different but dependent algebraic factors, which are de- 
termined by an electronic and nuclear coupling scheme in 
the upper hyperfine states and by the polarization of the 
exciting field. As an external characteristic of the cou- 
pling strength with the control mode, we use the Rabi 
frequency Oc = 2\Vnm'\/'h defined with respect to the 
lower energy transition. The results of our calculations 
will be compared with the usual A-scheme where only 
state \n) is taken into consideration. 

Figure [3] gives the spectral dependencies for imaginary 
(absorption) and real (dispersion) components of the sus- 
ceptibility x(^) (scaled in units of npA^) for a Rabi fre- 
quency J7c = 157 and for the control field on resonance 
with the atomic transition A = 0. In these conditions, 
there is no light absorption in the medium. The imagi- 
nary part of the susceptibility is responsible for the losses 
caused by incoherent scattering. All three AT compo- 
nents are shown as well as the doublet approximation 
calculated in the three-level A model when the upper 
state In') is disregarded. At first sight, the discrepancy 
between the exact result and the A-scheme approxima- 
tion seems small. However, two important qualitative 
differences can be pointed out. 

First, in the multilevel configuration, the spectral 
point, where due to electromagnetically- induced trans- 
parency (KIT) the absorption falls down to its minimum 
level, is shifted to the red and the KIT resonance does 
not occur for A = A, as it always does in the case of 
the A model. As clearly shown in the inset of figure [3) 
the light shift is significant and it depends on the Rabi 
frequency of the control field. This light shift was first 
observed by L. Hau and coworkers in the pioneering ex- 
periment on slow light [2T1. Second, the medium does not 
have a perfect transparency at the optimal EIT point as 
it also always does for the A-scheme approximation. In 
our example this effect seems rather weak but it is non- 
vanishing even in the limit of f^c ^ 0, and it can become 
important if the medium has a very large optical depth, 
where the losses caused by incoherent scattering can ac- 
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FIG. 3: (Color online) Absorption part (upper panel) and 
dispersion part (lower panel) of the sample susceptibility in 
the presence of the control field tuned in resonance with the 
atomic transition (A = 0), for a Rabi frequency Q.c = I57. 
The susceptibility components are scaled by noA^. Red solid 
curves show the result of the exact calculations and blue 
dashed curves are plotted for the usual A-scheme approxima- 
tion. For reference, grey curves give the profiles of the atomic 
resonances without control field. As shown in the inset of the 
upper panel, the EIT resonance is shifted to the red and the 
medium looses its perfect transparency at this point. 



cumulate [22' . 

In figure [4j the absorption and dispersion parts of the 
sample susceptibility are shown for two detunings of the 
control mode: A = — 5O7 and A — 5O7. For both de- 
tunings, the lines close to the atomic resonances are only 
slightly disturbed by the control field and we show only 
the part of the spectral domain where A ^ A. We com- 
pare the profile of the AT resonance located near the 
frequency of the control mode with predictions of the 
A-scheme approximation. The dependencies of figure HI 
show much larger deviation with the calculations based 



on the A-scheme than in the situation of figure |3] This 
is a direct consequence of the fact that the detuning A is 
comparable with the hyperfine interaction between elec- 
tronic angular momentum and nuclear spin. 

If the hyperfine interaction was negligible, which 
asymptotically occurs when one moves the AT resonance 
either much lower than state \n) or much higher than 
state |n'), then the considered scheme would be equiva- 
lent to the excitation of an atom with the nuclear spin 
decoupled from the electronic angular momentum, which 
has one-half value. The considered A-type interaction 
between the ground state Zeeman sublevels is impossible 
in this case. This explains why, when the detuning A 
is not negligible with respect to the hyperfine splitting 
^HF (here A = — 5O7 = — ri//i?/5.12), the amplitude of 
the AT resonance, shown in the upper panel of figure |4] is 
reduced as compared to the value given by calculations 
based on the A-scheme model. On the other hand, both 
models show approximately the same narrowing of the 
resonance line- width with increasing of A. More gener- 
ally modifications of the AT resonance in the vicinity of 
either Di or D2 transition due to the interplay between 
the hyperfine interaction and the optical detunings have 
a strong impact on observation of the EIT effect in the 
hot atomic gas and will be discussed in details in [52]. 

The situation changes if the control mode is tuned be- 
tween the hyperfine components as shown in the lower 
panel of figure |4] for A = +5O7. One can observe an 
enhancement of the AT effect in comparison with predic- 
tions of the A-scheme approximation. This enhancement 
is even slightly larger for the dispersion component, as it 
is more sensitive to the interference terms contributing 
into the sample susceptibility (see Eqs.(2.5) and (A4)). 



As a consequence the signal pulse passing through the 
atomic ensemble can undergo longer delay at the out- 
put than it would be expected from the A-scheme model. 
This constitutes an interesting advantage when the Ra- 
man process is applied to a quantum memory protocol in 
this configuration. This is detailed in the next section. 



IV. APPLICATION TO QUANTUM MEMORY 

As shown before, the Raman configuration with posi- 
tive detuning may lead to an improvement for the storage 
of a signal pulse. We evaluate here the performances of a 
memory in this particular case, where the control mode 
is tuned between the upper state hyperfine levels. In all 
the following, the Rabi frequency of the control mode is 
Q.C = 157 and it is detuned by A = 5O7, as in the lower 
panel of figure |4J 



A. Delay of the pulse transport in the atomic 
sample 

We consider a coherent pulse with rectangular profile 
impinging on the sample. Using Eqs.(2.5) and (2.31 we 
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FIG. 5: (Color online) Time dependence of the outgoing sig- 
nal pulses. The time i = corresponds to the arrival of the 
incoming pulse front to the atomic medium. For the reference, 
the gray box gives the outgoing pulse in the absence of atomic 
medium. Three pulses the with same duration T = 107~^ 
but different carrier frequencies a;_, uu and uu+ where ijj± are 
shifted by ±2n/T, propagate through the optically thick sam- 
ple with an optical depth uqX^L = 50. The inset shows the 
pulse spectrums with respect to the AT resonance. The tj_ 
pulse (red dotted line) experiences the longest delay but also 
the largest level of losses due to incoherent scattering. The 
tjj+ pulse (green dashed line) has the lowest level of losses but 
also the smallest delay. The central pulse has a intermediate 
carrier frequency and corresponds to a balance between the 
losses and delay. 



FIG. 4: (Color online) Same as in figurelslbut for two different 
detunings. The upper panel corresponds to A = — 5O7 and 
the lower panel to A = 5O7 . Only one component of the AT 
triplet, which is located near A ~ A, is shown. 



can calculate the shape of the outgoing pulse for var- 
ious situations. We have considered here three signal 
pulses having identical shape but different carrier fre- 
quencies. This approach allows to model a multimode 
situation where we make a discrete Fourier expansion 
with frequencies that are multiples of 27r/r, where T is 
the pulse duration. In the following, we use a dimension- 
less amplitude for the signal field denoted by a{z, t). We 
assume input pulses with rectangular profile, such that 
<^m{t) — 0{t) ~ 9{t — T), where 9{t) is the step- function, 
e{t) = l{t > 0) or = 0{t < 0). Figure [5] shows how 
this pulse is delayed and how its shape is modified after 
passing through the atomic medium. 

As shown in the inset of figure [5] the central pulse 
has a carrier frequency Hi chosen to ensure some bal- 
ance between transparency and delay. The two other 
pulses are shifted by ±2Tr/T, with Co± = a)±27r/T, which 
makes them spectrally orthogonal to the central mode. 
By extrapolating these series to infinite number of modes 



(w. 



2Trq/T, where q — 0, ±1, ±2, . . .) we can ex- 



tract the complete set of Fourier modes and use them for 
the expansion of an input signal pulse of arbitrary tem- 
poral profile with sharp bound. Such a Fourier descrip- 
tion of an arbitrary pulse of a finite duration T allows 
to follow the transformation of each spectral component 
in the output and could be important for the multimode 
quantum information processing. 

The delay of the pulse propagating through the op- 
tically dense medium of length L strongly depends on 
the optical depth. For a monochromatic signal mode at 
frequency ui the optical depth b(uj) is given by 



biCu) = A7rx"iCu) 



(4.1) 



The depth varies from b(uj) 3> 1 near the points of the AT 
resonances to b{ui) <C 1 in the transparency domains. In 
the A-scheme approximation, the resonance optical depth 
is defined in a unique way. In the multilevel situation 
the optical depth at resonance depends on the chosen 
transition. However it is always closes to bo ~ uqX L, 
which was taken equal to 50 for the round of calculations 
presented in figure |5] 

It can be seen in Figure [5] that the shape of the out- 
going pulses is strongly modified as compared to the ini- 



tial pulse. First, all three pulses corresponding to d), a)+ 
and W- are delayed. The pulse with frequency w_, the 
spectrum of which overlaps the most with the AT peak, 
spreads over longer times and has a much longer average 
delay than the other ones. However, this pulse has the 
highest level of losses due to the incoherent scattering. 
The pulse with carrier frequency aj+ has a lower level of 
losses and better preserves the original rectangular shape, 
but the delay is quite small. 

We can interpret the outgoing signal pulses in Figure 
[5] in the following way : the part of the pulse going out 
of the medium at times smaller than the time duration 
of the original pulse, T = 107^^, corresponds to direct 
transmission, while the tail of the outgoing pulse, which 
is emitted after the end of the incoming pulse, corre- 
sponds to a potentially stored and recovered signal. This 
is a good approximation if the transient processes asso- 
ciated with switching off/on the control mode as well as 
the retardation effects are negligible. The fact that in 
off-resonant stimulated Raman scattering the transient 
processes only weakly interfere with the transport of the 
signal pulse was verified in jTS] . 



B. Efficiency of pulse storage and retrieval 

We now turn to the full memory protocol, including 
storage and retrieval stages. The control pulse is switched 
off when the end front of the incoming signal pulse has 
entered the sample. After short transient dynamics, the 
signal wave packet, which is localized in the sample, is 
mostly converted into a standing spin wave distributed 
into the whole sample. This spin-wave is given by the off- 
diagonal matrix elements (spin coherence) between Zee- 
man sublevels m and m' existing in the system at time 
T, which we denote as a{z,T). To calculate the spin dis- 
tribution we have applied the theoretical approach de- 
veloped in [16' and modified it for the present case by 
including the dissipation processes. The spin distribu- 
tion is shown in the upper panel of figure [61 

After a controllable delay, the spin wave can be con- 
verted back into the signal wave packet by sending the 
second control pulse. This can be done in either forward 
or backward directions. In the first case, the retrieved 
pulse would reproduce the tail of the outgoing signal 
pulse, the propagation of which was interrupted in the 
writing stage of the memory protocol. This can be seen 
by comparing the graphs of figures [6] and [5] 

The pulse retrieval in the backward direction leads to 
some advantages. As discussed in [13 , while applying 
the time-reversal arguments one can expect a higher ef- 
ficiency for backward retrieval after a round of optimiza- 
tion. Also, the shape of the pulse recovered in the back- 
ward direction better reproduces its original profile. Our 
numerical simulations confirm this statement. The spin 
distribution for the mode uj- indicates that the signal is 
mostly stored near the left bound of the sample. Then 
it seems more natural to redirect the retrieved pulse in 
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FIG. 6: (Color online) Quantum memory. This figure shows 
the dynamics of the system for the same parameters as de- 
scribed in figure |5] Upper panel shows the spin distribution 
o{z, T) in the sample at time T ("stored light"). In the lower 
panel the signal pulse is retrieved on demand after a certain 
memory time via switching on the control field in either for- 
ward (solid lines) or backward (yellow filled area) direction. 
The leakage of each pulse, which is the light transmitted dur- 
ing the write-in stage, is indicated by the dotted line ((a): 
the input pulse carrier frequency is a;_, (b): the input pulse 
carrier frequency is uo ,(c): the input pulse carrier frequency 
is UJ+). The read-out stage is characterized by the memory 
efiiciency for the forward (fw) or backward (bw) retrieval. 



the backward direction where it would have less absorp- 
tion and, according to time-reversal arguments, would be 
regenerated by the medium to its original profile. That 
can be seen in figure [61 which shows approximately rect- 
angular retrieved profile for the a;_ mode with efficiency 



more than 60%. 

To conclude this section let us underscore the fact that 
the predicted quality of the considered memory channel 
is quite good. The obtained results show that in condi- 
tions currently attainable for a MOT the efficiency of the 
memory protocol could be expected up to 60%. We can 
expect further improvement in the efficiency by using the 
optimization of the temporal profile of the input pulse, 
as described in ref . [131 EH . 



V. CONCLUSION 

In this paper we have considered the propagation of 
a signal pulse in conditions of stimulated Raman pro- 
cess through a sample of alkali atoms. We have shown 
that the presence of the hyperfine interaction in the up- 
per state significantly modifies the Autler-Townes effect 
observed in the Di-line of such atoms. In particular, 
we have demonstrated that it is more efhcient to tune 
the control mode between the upper state hyperfine sub- 
levels where the stimulated Raman process is enhanced 
and provides more effective delay of the signal pulse. 

We have analyzed the dynamics of a signal pulse with 
initial rectangular profile and calculated how such a pulse 
is delayed and how its shape is modified after passing 
through the atomic medium. This dynamics reveals that 
the outgoing pulse is quite sensitive to the detuning of 
the pulse carrier frequency from the Autler Townes res- 
onance created by the control field. We have used our 
results to evaluate the efficiency of a quantum memory 
protocol based on the same stimulated Raman process. 
This has been done both for the forward and the back- 
ward retrieval schemes. We have shown that the back- 
ward retrieval can be more effective and the shape of 
recovered pulse can better reproduce its original profile. 
In particular, our results show that in conditions cur- 
rently attainable for a Cesium magneto-optical trap the 
efficiency of the memory protocol could be expected more 
than 60%. 



Maxwell equation considered in the main text of the pa- 
per can be expressed by the following diagram equation: 



= *^\J^\^ + 



*~\y\y^'\J\^ Jf\J^\f 



m 



(Al) 
Here the double wavy single ended lines describe the sig- 
nal field amplitude dressed by the interaction process 
shown in figure [l] when this field propagates through the 
sample. The single double ended wavy line is the retarded 
Green's function for light propagating freely in vacuum. 
The loop consisted of the atomic Green's functions and 
the vertices describe the polarization operator or suscep- 
tibility of the sample in response to the probe field. For 
the sake of clarity, we show the internal indices ni and 712 
running the quantum numbers of excited atomic states 
n, n' (figure [1]). 

The retarded-type Green's functions of the excited 
atomic states dressed by the interaction with the con- 
trol and vacuum modes and shown as double line in the 



polarization operator of Eq.(Al) are expressed by the 
following Dyson equations: 



+ - 




(A2) 
The outward and inward directed dashed arrows repre- 
sent the field amplitude of the control mode and the thick 
solid line is the retarded-type Green's function of state 
\ni) dressed by interaction with the vacuum modes only. 
This function in the first term in the right side should 
be multiplied by the factor (5„j ^^ (not shown) . It is ex- 
pressed by the diagram equation defining the natural de- 
cay rate of the excited state: 
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Appendix A: The atomic Green's functions and 
sample susceptibility 

The atomic Green's functions and sample susceptibil- 
ity are found in the second quantization formalism us- 
ing a non-stationary diagram technic [T51 [Ml [2S] . The 



-e ^ 



+ 



ii^l 



(A3) 
The thin solid and wavy lines are respectively non- 
perturbed vacuum Green's functions of the atom and 
field. 

The minus signs in these diagrams indicate the causal 
(time-ordered) character in averaging of the operators' 
product. In the case of field the expectation value is 
taken over its vacuum state. But in the case of atoms the 
averaging is taken over the initial state of atomic system. 
As a consequence, if in the above diagrams the time nat- 
urally rises from right to left, the backward directed thin 
line in the polarization operator ( Al I is proportional to 
the density of atoms in the initially populated state \m). 



The diagrams ( Al |-( A3 1 can be interpreted as a graphical 



solution of the Bloch equations in the first order response 
to the signal field. 

If the control mode is monochromatic then after 



Fourier transform the integral equations ( A2 ) can be 



written in a form of algebraic equations. These equations 
give us the set of four Green's functions Gnin2iP,E) de- 
fined in the reciprocal space. These functions contribute 



into equation (Al) with the "on shell" energy argument, 



which is given hy E = hioj + fi) + p^ /2mQ i. e. is ex- 
pressed by the kinetic energy of atom p^ /2mQ (p and toq 
are atomic momentum and mass respectively; the inter- 
nal ground state energy is considered as zero level) plus 
the energy of a photon from the signal field fi{u) + Vt). 
Let us point out that because of the momentum con- 
servation in the polarization operator the momentum of 
excited atom p' is shifted from p by the momentum of 
absorbed signal photon such that p^ = p^ -I- hk. 

The Green's functions have the following analytical ex- 
pression 



G^-^-\p,E) ^ hlE 
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Two other functions can be similarly written via trans- 
posing the indices n -s^ v! . 

The following notation is used: Vnm' = (dE)„^,, 
Vnira' — {'^^n'm' ^^^ '^'^ matrix elements of interac- 
tion with the control field; E is the complex amplitude 
of the field's positive frequency component E'^+)(r,i) = 
Eexp(— iwi 4- ikr); i?„, En' are non-perturbed energies 
of the excited states; 7 is the natural decay rate of the 
excited state. The quasi-energies in the denominators of 
Eqs.(A4) contribute to the excitation spectrum dressed 



by interaction with vacuum modes and control field while 
the hyperfine interaction is infinitely strong. They are 
given by 
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(A5) 



and similar for n' . Here Em' is the energy of state 
|m'), which for the system of degenerate Zeeman sub- 
levels shown in figure [l] coincides with the energy of state 
|m), such that Em' — Em = 0. Quasi-energies (A5) are 



convenient intermediate parameters describing the AT 
structure if only one state either \n) or \n') is taken into 
consideration. But actual location of the AT poles is 
described by full susceptibility of the sample given by 
Eq.(2.5), where all four Green's functions equally con- 



(A4) tribute. 
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